Guiding and Trapping of Electron Spin Waves in Atomic Hydrogen Gas 
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We present a high magnetic field study of electron spin waves in atomic hydrogen gas compressed 
to high densities of ~ 10 18 cm -3 at temperatures ranging from 0.26 to 0.6 K. We observed a variety 
of spin wave modes caused by the identical spin rotation effect with strong dependence on the spatial 
profile of the polarizing magnetic field. We demonstrate confinement of these modes in regions of 
strong magnetic field and manipulate their spatial distribution by changing the position of the field 
maximum. 

PACS numbers: 67.63.Gh, 67.30. hj, 32.30.Dx, 32.70.Jz 
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Spin and magnetization oscillations are well known 
phenomena in condensed matter physics. In solids wave- 
like spin excitations are caused by long-range dipolar 
or exchange interactions Quantized spin excitations 
(magnons) may exhibit collective quantum phenomena 
similar to BEC and superfluidity in quantum gases and 
liquid He. Observation of these effects in supernuid 3 He 
3 and ferromagnets Q have been extensively discussed 
recently. In gases spin waves may also propagate al- 
though interactions are weak and occur during short coi- 
lisional events. Quantum gases represent a special case 
where de Broglie wavelength A t h exceeds the character- 
istic range a of the interatomic potential and therefore 
the identical-particle exchange effects in atomic collisions 
appear long before the electrons of the atoms start to 
interact. If two colliding atoms have different spin direc- 
tions, exchange interaction leads to the rotation of both 
spins around their sum, a phenomenon called the Identi- 
cal Spin Rotation (ISR) effect Being accumulated 
in numerous collisions ISR leads to the propagation of 
spin excitations. The efficiency of ISR is determined by 
the parameter fj, ~ A t /j/et, which is large in the quantum 
gas limit. Therefore, for cold gases there exists a range 
of densities n and temperatures where n -1 / 3 3> A t h 3> a, 
i.e. the gas is in the quantum regime but not yet degener- 
ate. This range is especially wide for smallest atoms like 
H (a w 0.07 nm). Experimentally the ISR effect was first 
observed in the nuclear spins of electron-spin polarized H 
(HI) [6J], then in 3 He gas [7j and recently also in ultracold 
87 Rb vapors [a. M ■ ISR waves were predicted for electron 
spins of HJ, [J, [l(| but have not been observed so far. In 
this letter we report a quantitative study of electron spin 
waves in HJ, and verify that they are associated with ISR. 

Spin transport in quantum gases is described by 
the complex diffusion equation for the transversal spin- 



polarization S-i 
and fi 3> 1 to: 
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where Dq is the spin diffusion coefficient in unpolarized 



gas, e = +1 for bosons and —1 for fermions, 7 is the 
gyromagnetic ratio, and SBq is the deviation of mag- 
netic field from its average value Bq. Eq. ([TJ is similar 
to the Schrodinger equation for a particle with an effec- 
tive mass M* = —h(i/2DoS and potential energy given 
by the second term on the right-hand side. The behav- 
ior of such quasiparticles in an inhomogeneous magnetic 
field depends on the signs of the parameters fj,, s, and 7. 
The spin-wave quality factor /i depends on the details of 
the interatomic potential. For H fi < 0, and for 3 He it 
changes sign from negative to positive when temperature 
is lowered below w 0.5 K From Eq. (JXJ) we see that 
similarly to low-field seeking atomic species the nuclear 
spin waves of H have lower potential energy in regions 
of smaller magnetic field. Due to 7 being negative the 
electron spin waves behave as high field seekers. Simi- 
lar behavior is expected for 3 He at T > 0.5 K (/i > 0, 
e = —1). In addition, as 7 e /7n ~ 650, the strength 
of the potential energy for electron spin waves is much 
larger than that for nuclei. Consequently local magnetic 
field maxima become effective potential wells. We will 
show in this letter that similarly to real particles ISR 
spin waves in HJ, may be spatially confined and manipu- 
lated by magnetic forces. 

In our experiment HJ, gas is hydraulically compressed 
up to n ~ 10 18 cm~ 3 , as described in The sample 
is located in a thin- walled (2 /im) 0.5 mm diam. Kapton 
tube (KT) below the flat mirror of a Fabry-Perot res- 
onator (FPR) (see Fig. [TJ. The sample volume (SV) is 
coupled to the rf. field of the resonator via an evanescent 
field (EF) beneath a subcritical (0.4 mm dia.) orifice in 
the FPR-mirror. The evanescent field amplitude Hi(z) 
has an approximately Gaussian shape decreasing down- 
wards, with an effective height of Ief ~ 80 /im [l4( . The 
bottom of SV is limited by the meniscus of supernuid 
helium. The height L of SV depends on the position of 
the helium level and can be varied from 6.5 to 0.5 mm. 
For smaller L the sample evolves into a bubble. The 
Kapton tube is glued with Stycast 1266 to an epoxy disk 
with conical cross section. The thickness of the disk at 
the tube wall is rj 0.5 mm. The outer surface of KT 
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FIG. 1: (Color online) (a) Schematic drawing of the sam- 
ple cell showing the sample volume (SV), Kapton tube (KT), 
liquid- helium volumes (1-He) and epoxy structures (E). (b) 
Magnetic field inhomogeneity due to magnetized epoxy. The 
field amplitude in the top of SV is shown as a colormap, black 
(red) corresponding to the strongest and blue (white) to the 
weakest field. Axial (z) and radial (r) cuts of the field profile 
through the field maximum are also shown for three differ- 
ent additional linear gradient values: -1.7 (red/black), -0.85 
(green/dark) and (blue/light) G/mm. 



below the epoxy disk is flushed with superfluid helium 
cooled to 200-600 mK in a separate heat exchanger cou- 
pled to the mixing chamber of the dilution refrigerator. 
This construction allows to reach high densities with lin- 
ear response of the FPR to the electron spin resonance 
(ESR) in HJ,. For the excitation and detection of spin 
waves a highly inhomogeneous rf. held combined with 
a controllable static field are used. The homogeneity of 
the magnetic field in the compression region was lim- 
ited to AttM ~ 0.8 G by the weakly paramagnetic epoxy 
disk, the field of which is equivalent to that of a minia- 
ture solenoid located at the tube wall, creating a saddle 
shaped field with a maximum located near the wall (see 
Fig. [lib)). We were able to reduce the inhomogeneity 
to ss 20 mG over the rf. field region using a set of linear 
and parabolic shim coils, which were also capable of gen- 
erating linear axial field gradients up to G z — ±4 G/mm. 

The H.J. samples were studied by ESR technique at 128 
GHz [l5| . The spectra were recorded by applying a con- 
stant frequency from a highly stabilized mm- wave source 
while sweeping the static magnetic field offset h through 
the resonance. We utilized the b-c transition, since the 
population of the hyperfine state a is vanishingly small 
due to rapid exchange recombination (for the notations 
of the hyperfine states see, e.g. [H}). At small enough 
densities n < 5 x 10 16 cm~ 3 the ESR lines were inho- 
mogencously broadened due to the spurious field of the 
epoxy ring. A small positive gradient G z ~ 1 G/mm best 
compensated the inhomogeneity resulting in the smallest 
line widths. For n > 10 17 cm~ 3 we observed two distinct 
changes in the shape of the ESR lines: (i) the main ab- 
sorption peak was split into several narrow lines whose 
position and separation depended on n and G z , and (ii) 
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FIG. 2: ESR spectra for different values of axial magnetic field 
gradient. The resonance frequency is fixed and magnetic field 
is swept through the resonance. The regions of the sample 
volume located in larger local fields will get in resonance at 
smaller sweep field. The direction of the field axis is reversed 
to get the local field increasing from left to right. Insert: peak 
(ii) position as a function of G z . 



within the range of G z s» —2... + 0.85 G/mm an extra 
peak appeared on the right side of the spectra. These 
features are summarized in Figs. [2] and [3l The separa- 
tion of this peak (ii) from the main ESR line increased 
with G z and remained visible even at the lowest densi- 
ties n ~ 10 16 cm -3 at which we were still able to resolve 
ESR lines. However, the peak (ii) was never observed on 
the left from the main resonance line not even with suf- 
ficiently large positive gradients. Neither was there any 
dependence of the peak positions on the height of the gas 
sample, except for L < 0.5 mm, when the sample evolved 
into a small bubble. 

Next we consider the splitting of the main ESR peak 
for high densities and large negative G z (feature (i)). To 
estimate the expected separation between the ISR spin 
wave modes we assume an integer number of their half- 
wavelengths between the walls of the compression region. 
From the ISR spin- wave dispersion law uj — ^jf-k 2 we get 
Alu < tt 2 Dq/^,L 2 . Taking the minimum sample height 
L w 0.5 mm and the values D^n = 1.5 x 10 18 cm _1 s _1 
and u k 7 from experiments with nuclear spin waves in 
HJ. [6j, we find that the spin wave modes will be sepa- 
rated by < 1.5 kHz at n — 10 17 cm~ 3 . This corresponds 
to < 0.5 mG in the field sweep and cannot be resolved in 
the experiment. It is thus unlikely that the peaks seen in 
Fig. [5] and [3] correspond to individual spin wave modes, 
but rather to a number of overlapping modes. It turns 
out that just in the case of strong negative field deviation 
SHq = G z z the overlapping modes provide a resolvable 
structure. The small inhomogeneity created by epoxy 
may be neglected in this case. A negative G z implies 
that the potential shape is a box with tilted bottom de- 
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FIG. 3: ESR spectra observed for samples G z = —4 G/mm. 
Dotted line: the simulated contribution F(h) of spin waves 
to the ESR absorption. In the insert the orientations of the 
potential U(z) (thick solid line), a spin wave mode S+(z), and 
the evanescent field Hi(z) are shown emphasizing the upper 
part of the sample volume. 

creasing downwards as U(z) = 7 e G z z (insert in Fig. [3]). 
The oscillating solutions of Eq. (p} are found in the form 
of Bessel functions in the radial and Airy functions in the 
axial direction. The shape of the solution should match 
the shape of the rf. field excitation, which is axially sym- 
metric and has a maximum at r = 0. Since only the 
zcroth-order Bessel function satisfies this condition, the 
solution is written as 

S+ ~ J (Kr)e iut Ai(z/\ + \ 2 k 2 z ), (2) 

with k 2 = (ilj/Do and A = (D /7eG z /u) 1//3 being the 
length scale of the Ai(z) function. We apply reflective 
boundary conditions with zero spin flux through the walls 
of the cylinder. These conditions define the eigenvalues 
of k and k z , eigenfrequencies oj of the oscillations, and 
the number of Ai function oscillations in the axial direc- 
tion. For n = 5x 10 17 cm~ 3 and T — 260 mK we esti- 
mate A w 7 /*m. In order to generate spin waves, their 
wave function given by Eq. ([2]) should overlap with the 
rf. field. This implies that in a negative gradient we may 
excite only high-lying modes, which have their frequency 
eigenvalues in the vicinity of the intersection of the linear 
potential and the upper wall (see insert in Fig. [3]). The 
total number of modes counting from the bottom of the 
potential U(z) and the corresponding number of Ai(z) os- 
cillations are 3> L/X ~ 100 even for our shortest samples. 
Therefore, we deal with a quasi-continuum of mode fre- 
quencies, and excite only those with non- vanishing wave 
functions near the top of the cylinder. The position of 
the lower boundary does not influence the shape of these 
modes near the top because of the rapid oscillation of 
Ai(z) near the bottom. This explains why the observed 
ESR lineshapes do not depend on the sample height. 
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FIG. 4: Positions of magnon peaks for samples in large neg- 
ative magnetic field gradients of -4 G/mm as functions of H 
gas density. 



Since the position of the lower boundary does not influ- 
ence the shape of 5-1- near the top, we consider the case of 
a semi-infinite cylinder. In this case the spin waves may 
freely propagate downwards along the slope of the mag- 
netic potential. Sweeping the magnetic field h we scan 
through the region Ief by tuning the rf. field in reso- 
nance with electron spins. At each point the transversal 
magnetization is created locally and two spin waves are 
launched; one up and the other down. For zero axial gra- 
dient these would be similar to the electromagnetic waves 
in a cylindrical waveguide. Moving the point of excita- 
tion with respect to the top we scan a standing wave 
pattern formed by the direct and reflected waves with 
the spatial period equal to half- wavelength. For non-zero 
gradients the sine-waves are accelerated downwards and 
are replaced by the Airy functions Ai(z/X). The contri- 
bution of spin wave modes to the ESR absorption is then 
~ F(h) = Ai 2 (/i/(G z A)) x H 2 (h/G z ), shown by the dot- 
ted line in Fig.[3j This will produce maxima of absorption 
at the antinodes of Ai(z/X), the points Zj corresponding 
to the roots Xi of dAi(x)/dx. The peak positions in the 
field sweep are hi = G z Xxi + Bq. They have a n -1 / 3 den- 
sity dependency since the diffusion coefficient is inversely 
proportional to the H density; nD = 1.5 x 10 18 cm -3 
[l2|. The separation between the peaks is proportional 
to the roots of the Ai function, which we verified as fol- 
lows. First we fit the position of the leftmost peak by the 
relation hi — Gn~ x ^x\ +Bo with C and Bo being fitting 
parameters. Then we calculate the positions of the next 
four peaks as functions of density as hi = Cn~ 1 / 3 Xi + Bo- 
The results plotted by dashed lines in Fig. @] match the 
experimental data within the error bars. The relation 
C = (G 2 nDo/ '7 e ^) 1 / 3 provides the measurement of the 
ratio nDo/fi = 1.9(4) x 10 17 cm -3 . From the fit we find 
H = 8(2), which has been confirmed by numerical solu- 
tions of Eq. ((T|) . The value of /i obtained in our work is 
in good agreement with the results for nuclear spin waves 
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in Hi |6j. 

The feature (ii) was observed for small or zero linear 
gradients. In this case there is a local magnetic field 
maximum near the top of the sample and close to EF. 
An example of a SB(r, z) profile calculated for G z = is 
presented in Fig.[T](b). Together with the wall of the tube 
such a field profile forms a 3D trap for the electron spin 
waves. A numerical solution of Eq. ([I]) reveals a large 
number of modes in this trap with the separation of a 
few kHz, again too small to be resolved. Sweeping the 
static field offset h we tune our excitation to resonance 
in different regions of the sample volume. Absorption is 
increased as energy is pumped into the spin wave modes 
and the absorption strength depends on the density of 
modes for a given frequency and spatial position. Ap- 
plying a small linear gradient we modify the axial profile 
and consequently move the trap bottom up or down as 
shown in Fig. Q] (b). The position of the corresponding 
spin wave peak in the ESR spectrum is determined by the 
static field at the trap bottom, which leads to a nonlin- 
ear dependence of the peak position on G z (Fig.[2]insert). 
This also explains why we cannot see the feature (ii) on 
the left side of the main ESR line. For large positive 
gradients the field maximum is always located at the top 
of the cylinder, coinciding with the position of the main 
ESR line at the H^(z) maximum. Note that if the peak 
(ii) would be caused by the second rf. field maximum lo- 
cated some what below the main one, changing the sign 
of the field gradient would mirror the spectrum and this 
peak (ii) would be on the left hand side with linear de- 
pendence on G z . Remarkably, the spin- wave modes (ii) 
are well excited and detected even by the very weak tail 
of the rf. field. The spin wave signal is much larger than 
absorption due to the ESR itself. Comparing the density 
dependence of the two spin wave we find that the spin 
waves (ii) survive at much smaller gas densities than the 
propagating modes (i). These observations confirm that 
in case (i) we deal with a spin wave-guide and with a 3D 
spin wave trap in case (ii). 

We have also considered the possibility of spin exci- 
tations called magnetostatic Walker modes [17| which 
should exhibit a strong dependence of the mode posi- 
tions on the sample geometry. However, this was not 
observed in our experiments. The magnetostatic modes 
do not depend on temperature, because they are caused 
by long range dipolar forces and do not involve atomic 
collisions. We verified that the separation between peaks 
(i) increases slightly at higher temperatures. For the ISR 
modes we expect a temperature dependence Do /m ~ VT 
or weaker [12} , which is in line with our observations. 
We conclude that both types of the spin wave modes de- 
scribed above are caused by the ISR effect. 

In this work we have studied electron spin waves in 
atomic hydrogen quantum gas. We detected two types 
of spin wave excitations; travelling modes guided by the 
cylindrical spin wave-guide and modes confined in the 



magnetic potential well. Similarly to the matter wave in- 
terference of cold atoms [HI, [l9[ the presence of a reflec- 
tive boundary in the spin- waveguide leads to interference 
of the two travelling spin waves, one falling straight down, 
and the other one reflecting from the upper wall. In the 
quantum regime the spin waves are described as quasi- 
particles called magnons 0, Q . Our trapping technique 
resembles the original technique for trapping cold gases, 
but here it is realized for high field seeking quasiparti- 
cles. In further experiments we hope to observe effects of 
statistical correlations between trapped magnons leading 
to Bose-Einstein condensation and spin superfluidity. 
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